Correlation driven Mott-insulator-to-metal transition 
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We study transport properties of the half-filled two-dimensional (2D) Hubbard model with spa- 
tially varying interactions, where a pattern of interacting and non-interacting sites is formed. We use 
Determinantal Quantum Monte Carlo method to calculate the double occupation, effective hopping 
and Drude weight. These data point to two phase transitions, driven by fermionic correlations. The 
first is the expected metal to a Mott insulating state. The second one, is an exotic transition from a 
Mott insulating state to a highly anisotropic metal, that takes place at large values of the fermion- 
fermion interaction. This second transition occurs when the layers formed by the spatially varying 
interactions decouple due to the suppression of the hopping between interacting and non-interacting 
sites, leading to fermionic transport along the non-interacting ones. 
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PACS numbers: 71.10.Fd 71.30.+h 71.27.+a 

Metal insulator transitions have been a topic of intense 
study over many years 0. When the transition from a 
metallic into an insulating state is driven by correlations 
it is known as the Mott transition 0, H| ■ The experimen- 
tal observation of the Mott insulating state in two-flavor 
mixtures of fermionic atoms loaded on optical lattices, 
0, HI opened new possibilities in this field. Models that 
take correlation effects in tight-binding systems, such as 
the Hubbard model, have been widely used to describe 
different types of insulating phases and metal-insulator 
transitions. 

An interesting point that has been addressed in the 
literature is whether these same correlations could drive 
an insulating system metallic. When a periodic on-site 
energy is distributed in the system, we have the so called 
Ionic Hubbard Model. In the non-interacting limit, the 
periodic potential produces a dispersion relation that is 
gapped at half-filling. In these systems, correlations are 
responsible for a band-insulator to metal transition. 

0- 

!■ 

The same on-site energies can be randomly distributed 
thr oug hout the lattice, leading to an Anderson insulator 
H> El depending on the strength of the disorder and 
system dimensionality. The use of speckle noise (Tl. 12\ 
in cold atoms loaded on optical lattices has enabled a 
new experimental route to understanding the interplay 
between disorder and interactions. This is still an area of 
intense debate, with many open questions. Nonetheless, 
numerical data for the Anderson-Hubbard model points 
to an Anderson insulator to metal transition driven by 
correlations 13 1. 



Can correlations drive a Mott insulator into a metal? 
Here we show that, when the interactions are not ho- 
mogeneously distributed in the system, but instead are 
spatially varying, this correlation induced Mott insulator 
to metal transition can indeed take place. Although spa- 
tially varying interactions are not yet available on optical 
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FIG. 1: (Color online) Schematic figure of the spatially vary- 
ing interactions in a system with Lu = 2 Lo — 1 in a 12 x 12 
lattice. 



lattices, they have recently become available in ultracold 
gases. The ability to control a magnetic Feshbach reso- 
nance with laser light [ljj], has increased the tunability 
of interactions for bosonic systems. Submicron spatial 
modulation of the interaction was already achieved in a 
174 Yb Bose-Einstein condensate (ljj . Optical control of 
Feshbach resonances for a fermionic ultracold gas [l6| has 
also been proposed. 

Here we consider a modified version of the Hubbard 
Model (HM) with spatially varying interactions, whose 
Hamiltonian reads 

n = -* E 

i ^ ' ^ ^ ' No- 

where, in standard notation c\ a (c ia ) creates (destroys) 
a fermion in site i in state a. The hopping parameter 
between nearest neighbor sites is set to t = 1, 

fTj is the site dependent interaction, and is the chemi- 
cal potential controlling the band filling to yield a given 
fermionic density p. The interaction term is written in 
particle-hole symmetric form. Thus, tuning /i = drives 
the occupation to one in every site for all Hamiltonian 
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FIG. 2: (Color online) Double occupation within (a) repulsive 
(du) and (b) non-interacting sites (do) as a function of U for 
different 12 x 12 superlattices at T = 0.10. Double occupa- 
tion for the corresponding homogeneous system (stars) is also 
shown in (a) for comparison. 



parameters t, Ui and temperatures T . We have here re- 
stricted our study to half-filled systems. To simulate the 
spatially modulated systems we construct a layered pat- 
tern of repulsive and non-interacting layers where U > 
and U = 0, respectively. We define the width of the re- 
pulsive layer as Ljj and that of the non-interacting one 
as Lq as depicted in Fig. [TJ for a Ljj = 2 Lq = 1 in a 
12 x 12 lattice. Note that not all patterns fit all the avail- 
able lattice sizes. Here we use Determinantal Quantum 
Monte Carlo (DQMC) simulations [lj| to probe trans- 
port properties of the half-filled two-dimensional square 
L x L system with spatially varying interactions. Metal- 
insulator transitions of one-dimensional non-symmetrical 
Hubbard superlattices [171 ] were studied and shown to 



have an interesting behavior where the Mott transition 
takes place away from half-filling, at a pattern dependent 
density 0. 

A quantity of interest is the double occupation d = 
(n-j-nj,), as it is measured in optical lattices experiments. 
Indeed, it has been used to probe metal insulator transi- 
tions [5| in fermionic two-flavor mixtures of 40 K atoms. 
For the homogeneous, non-interacting, metallic system 
d = (nf)( n 4-) = 0.25. As correlations are turned on at 
a fixed temperature, d decreases, signaling the metal- 
insulator transition, and approaches as U — >• oo. The 
symmetric form of the Hamiltonian requires that, at half- 
filling, the charge distribution is homogeneous through- 
out the lattice. The double occupation, on the other 
hand, is not homogeneous and follows the superlattice 
pattern. Single-atom resolved double occupation [20] is 
experimentally available, therefore it is then relevant to 
define do an d djj, the average double occupation within 
non-interacting and interacting sites, respectively. 



The [/-dependence of both do and djj for different su- 
perlattice patterns and fixed temperature (T = 0.10) is 
shown in Fig. [5J Within the repulsive sites (Fig. Ufa)) 
du decreases monotonically with U, and approaches that 
of the corresponding homogeneous system as the ratio 
Ljj I Lq is increased. The behavior of djj is in agreement 
with a metal-insulator transition as the interactions are 
turned on. The double occupation on the non-interacting 
sites (Fig. H{b)) is affected by the strength of the in- 
teraction on the neighboring repulsive ones and displays 
a surprising non-monotonic behavior with U. Starting 
from U = 0, as U increases the effect of the interactions 
"leaks" into the non-interacting sites. When the ratio 
Ljj /Lq increases, the double occupation is reduced even 
where the interaction is turned off. It is worth noting 
the difference in the cases Ljj = 1 Lq = 1 and Ljj = 2 
Lq = 2, where the ratio Ljj/Lq is the same, but the wider 
free layer in the latter makes the spins more uncorrelated 
in the non-interacting region. As U — » oo double occupa- 
tions are forbidden in the repulsive sites and hopping be- 
tween free and repulsive sites is complctclly suppressed, 
leading to a decoupling between free and repulsive layers 
and bringing do back to the value of the homogeneous 
non-interacting system, do — > 0.25. The upturn in do 
signals the onset of the decoupling between the layers, 
that, as we shall see bellow, leads to an insulator-to-metal 
transition. 

, which 



We now analyze the effective hopping [21 



is defined as the ratio of the kinetic energy on a super- 
lattice to its non-interacting homogeneous counterpart 
value, both in the direction along (x) and across (y) the 
layers; 



t, 
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(2) 



where a = x or y. One would expect that anisotropy 
favors the fermionic transport along the direction of the 
layers. This is indeed the case, as can be readily observed 
in Fig. [3J where the repulsive interaction splits the x 
and y contributions to the effective hopping, reducing 
ty.eff/ty more strongly than t Xte ff/t x as U is increased. 
In the direction perpendicular to the layers the effective 
hopping approaches that of the corresponding homoge- 
neous system, as Ljj/Lq increases, displaying the char- 
acteristic behavior of an insulating system. Along the 
direction of the layers, on the other hand, the effective 
hopping is close to that of the non-interacting metallic 
system when Ljj/Lq < 1. 

We can gain further insight by analyzing the effective 
hopping separately within free and repulsive layers. This 
is done in Fig. |H where t^ e j,f/t x (filled symbols) and 
t° x e ff/t x (empty symbols) are normalized by the number 
of repulsive and free sites, respectively. The effective hop- 
ping along the repulsive sites is independent of the super- 
lattice pattern for small U and follows closely the analyt- 
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FIG. 3: (Color online) Effective hopping along (a) and across 
(b) the direction of the layers as a function of interaction 
strength U, for 12 x 12 superlattices at T — 0.10. Results for 
the corresponding homogeneous system (stars) are shown in 
both panels for comparison. 
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FIG. 4: (Color online) Effective hopping along free (open sym- 
bols) and repulsive (closed symbols) for 12 x 12 superlattices 
at T = 0.10 as a function of U. Dashed line is perturbation 
theory analytical result, from Ref. [2l|. 



ical results for perturbation theory in the homogeneous 
system [2lj ■ As U is increased there are small deviations 
from the homogeneous counterpart. From the similar- 
ity of the effective hopping in the homogeneous system 
and along the repulsive layers, we can infer that trans- 
port properties are essentially the same in both, leading 
to an insulating behavior of the fermions along the cor- 
related layer. Therefore, the main contribution to the 
transport in the direction parallel to the layers comes 
from the non-interacting sites, where t^. e f*/t x > 1. Fig. 
5] shows that as U is increased the contribution to the 
kinetic energy due to hopping between free sites is en- 
hanced. When comparing different patterns we can see 
that the enhancement in t° x e ^ is optimized in lattices 
with narrow non-interacting layers (Lq = 1) separated 
by wide repulsive ones. The picture emerging from Fig. 
|4]is that, while the repulsive layers are insulating, trans- 
port along the non-interacting layers is actually improved 
by the presence of the neighboring repulsive sites and in- 
creasing interaction strength among them. 

The effective hopping cannot be used as an order pa- 



rameter for the metal-insulator transition as it only ap- 
proaches zero for U — > oo in the homogeneous case, 
whereas the system is known to be in a Mott insulating 
state for any non-zero value of U. Among the available 
quantities to characterize the transport properties of the 
HM within DQMC simulations, the Prude weight is the 
one less affected by finite size effects 
termincd by 24 1 



23j. It can be de- 



D P . 



lim \{-k a ) - A QQ (q = 0, iu m )], 

T.m—^Q 



(3) 



where (k a ) is the contribution to the average kinetic en- 
ergy from fermions hopping along the a direction and 
A QQ (q = 0, iuj m ) is the long wavelength limit of the 
current-current correlation function with u m = IrrnxT 
the Matsubara frequency and a, = x or y. In principle 
the Drude weight as defined by (j3]) only provides infor- 
mation about the ground state of the system. Here we 
calculate temperature dependent approximants to D and 
analyze how they converge as the temperature is lowered. 

The Drude weight takes different values along the two 
preferential directions. Across the layers it displays an 
insulating character for all superlattices considered (not 
shown) in agreement with the effective hopping behavior. 
In Fig. [5ja) we plot the Drude weight along the direction 
of the layers for a fixed system size (8 x 8) for different 
patterns as wells as for the homogeneous system. The 
Drude weight for the U ^ homogeneous system (stars) 
shows the expected insulating behavior; as the tempera- 
ture is reduced D/ire 2 — > at a given temperature, we 
shall call T*. For U = 4 (dot centered stars), T - 0.2 and 
for U = 10 (closed stars), T* - 0.85. T* is related to the 
charge gap and increases as U increases. Our values of 



T* are in agreement with the results of Ref. |25j where this 
temperature is obtained in a completely different fashion, 
by the gap emergence in the spectral weight function. 
The negative values of D are a feature of the finite sys- 
tem size 24j , and converge to zero in the thermodynamic 
limit [23 1 . For the non-interacting homogeneous system 
(open squares) D/ire 2 = 0.79 as T — > 0, indicating an 
isotropic metallic behavior. 

Let us now turn to the superlattices, with an interme- 
diate interaction strength U = 4 (dot centered symbols), 
the Drude weight goes to zero, signaling an insulating 
ground state; T* depends on the interaction pattern and 
increases as Ljj/Lq increases, being smaller than its ho- 
mogeneous counterpart. Contrary to what happens in 
the homogeneous system, increasing U does not increase 
the charge gap, see closed symbols in Fig. 03a). For the 
system with Ljj = 3 and Lq = 1, T* is roughly the same 
for U = 4 and U =10, whereas for the systems with 
Ljj/Lq < 1 the Drude weight stabilizes at non-zero val- 
ues for the lower temperatures achieved, and approaches 
that of the non- interacting case as L\j / L is reduced. 

In Fig. [S^b) we study size and temperature effects 
on the Drude weight for the Ljj = 1 Lq = 3 system. 
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FIG. 5: (Color online) (a) Temperature dependence of the 
Drude weight along the direction of the layers, for 8x8 lattices 
for different SL's with U = 4 and U = 10 and homogeneous 
system with U = 0, 4 and 10; (b) U dependence of the Drde 
weight for the Ljj — 1 Lq = 1 superlattice at T = 0.0625 and 
T = 0.1000. 



As D x /ire 2 evolves with U we can divide the system's 
behavior in three regions, determined by the crossings 
between 8x8 and 16 x 16 data. Let us call the val- 
ues of U where the crossings take place Uci and Uci- 
For U < Uci one can see that D x /ire 2 is larger for the 
larger system size, pointing to a non-zero contribution as 
T — > and L — > oo, the characteristic metallic behavior. 
For Uci < U < Uc2 we have the opposite behavior; as 
the system size is increased, the Drude weight decreases, 
leading to a vanishing contribution in the thermodynamic 
limit. Note that this insulating region increases as the 
temperature decreases. As U increases even further and 
the region beyond Uc2 is reached, the behavior is similar 
to that of the first region: D x /ire 2 is larger for the larger 
system size, and the system is once again metallic. The 
values of Uci and Uc2 depend on temperature and are a 
likely to change when larger system sizes are considered. 
Nonetheless, the existence of two critical U values lead- 
ing to a metal to Mott insulator transition for small U, 
followed by a transition from an insulating to a metallic 
state along the direction of the layers at large U is sup- 
ported not only by the Drude weight data, but also from 
the effective hopping and double occupation. 

In summary, we have employed Determinantal Quan- 
tum Monte Carlo simulations to study the effect of 
spatially varying interactions in two dimensional half- 
filled systems strongly interacting fermions. We have 
found that, as expected, when interactions are turned 
on the system undergoes a metal to Mott insulator tran- 
sition. Although the system is particle-hole symmetric 
and charge is evenly distributed among free and repulsive 
sites, the double occupation is not constant throughout 
the lattice. Double occupation on non-interacting sites 
are influenced by the neighboring repulsive sites, first de- 
creasing as U increases and then increasing to return to 



the uncorrelated value as U — > oo. This upturn in the 
double occupation signals the decoupling between free 
and repulsive layers that takes place at large values of U. 
Evidence from double occupation, effective hopping and 
Drude weight show that, when this decoupling occurs, 
the system undergoes an exotic transition from a Mott 
insulator to metal, driven by fermionic correlations. We 
hope this results will encourage the experimental real- 
ization of such spatially varying interactions in optical 
lattices. 
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